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Abstract. The hydrodynamic response of the inviscid small shearing box model of a midplane section of a 
rotationally supported astrophysical disk is examined. An energy functional £ is formulated for the general 
nonlinear problem. It is found that the fate of disturbances is related to the conservation of this quantity which, 
in turn, depends on the boundary conditions utilized: £ is conserved for channel boundary conditions while it 
is not conserved in general for shearing box conditions. Linearized disturbances subject to channel boundary 
conditions have normal-modes described by Bessel Functions and are qualitatively governed by a quantity E 
which is a measure of the ratio between the azimuthal and vertical wavelengths. Inertial oscillations ensue if 
E > 1 - otherwise disturbances must in general be treated as an initial value problem. We reflect upon these 
results and offer a speculation. 
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"Extraordinary how mathematics help you to know 

yourself." 
Samuel Beckett, Molloy 

1. Introduction 

The hydrodynamic response in small shearing box (SSB 
hereafter) models of rotationally supported flows has 
been receiving renewed attention over the past few years. 
Because such fluid configurations do not support su- 
percritical linear instabilities alternative forms of transi- 
tion, possibly subcritical, have been proposed, discussed 
and explored in numerous recent studies (Ioannou & 
Kakouris, 2001, Longaretti, 2002, Chagelishvilli et al, 
2003, Yecko, 2004, Mukhopadhyay et al 2004, Afshordi et 
al. 2004, Umurhan & Regev, 2004, Sternberg, 2005, Lesur 
& Longaretti, 2005, Balbus & Hawley, 2006, Rincon et 
al. 2007, Ji et al, 2006, Williams, 20063, Jimenez et al, 
2006, Shen et al, 2006). The identification of a route to 
sustained nonlinear activity for three dimensional distur- 
bances has been elusive for flows that are wall-bounded 
(see recently Rincon et al 2007) and for flows having open 
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boundaries (Shen et al, 2006) although weak levels of ac- 
tivity have been reported by Lesur & Longaretti (2005) 
near the Rayleigh line in flows with open boundaries. 

Far from answering whether or not such transitions do 
exist in these flows, we instead identify two clues which 
could be used to cast some light onto this question in the 
future. To this end we reconsider the dynamics of invis- 
cid anticylonic Rayleigh stable flows in the SSB model 
(i.e., rotating plane Couette flow, rpCf hereafter) subject 
to shearing boundary conditions (Goldreich & Lynden- 
Bell, 1965, Rogallo, 1981, Knobloch, 1984, Korycansky, 
1992, SBC hereafter) or channel boundary conditions 
(e.g., Yecko, 2004, Mukhopadhyay, et al. 2004, Lesur & 
Longaretti, 2005, Rincon et al. 2007, CBC hereafter). 

The first result is a general nonlinear feature pertain- 
ing to the energetics of such flows. If one defines £ as the 
total energy of the flow in the SSB, whose dynamical con- 
stituent includes the total kinetic energy, (i.e., the energy 
formed by considering the total velocity - which is to say 
the velocity fluctuations plus background shear flow) , the 
fate of £ depends upon the boundary conditions employed. 
£ is conserved for (but not limited to) (i) disturbances 
subject to CBC, (ii) azimuthally and vertically periodic, 
radially localized disturbances in radially unbounded do- 
mains. On the other hand, disturbances subject to SBC 
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do not identically conserve £ and its temporal quality de- 
pends critically upon the dynamics within the domain. 
The nonconservation of £ for flows subject to SBC relates 
to the feature that the kinetic energy of a fluid parcel 
jumps when it crosses the (quasi) periodic radial bound- 
aries. 

For our second result, we find that there are qualita- 
tive differences in the behavior of linearized disturbances 
subject to CBC depending on whether or not the quantity 
E is less than or greater than 1. This parameter is defined 
as 



where £ z ,£ y are the disturbance length scales in the ver- 
tical (like the disk normal) and streamwise (like the disk 
azimuth) directions of any given modal disturbance. The 
global flow for which the SSB is meant to represent is char- 
acterized by a rotation rate Q that depends on the radial 
coordinate R according to the relationship Q, cx R~ q . This 
means that on the local scale there is a Couette veloc- 
ity profile given as U = —qflaxy. Keplerian flow is when 
q = 3/2. Additionally and hereafter, we will consider lin- 
ear shear profiles that are anti-cyclonic, i.e., those in which 
q > 0. At infinite Reynolds numbers (Re = oo) axisym- 
metric linear instability (Taylor vortices) occurs for q > 2. 
The significance of £ was recognized in previous studies 
(Knobloch, 1984, Korycansky, 1992, Sternberg, 2005, as 
the parameter (3 in Balbus & Hawley, 2006) of distur- 
bances subject to SBC. In those studies it was demon- 
strated that modes always decay. If the mode was initially 
leading then it will experience one maximum whereas if 
the mode is trailing it never experiences a maximum. 
However, irrespective of this the asymptotic decay of dis- 
turbances depends on S: for £ < 1 the decay is monotonic 
while for S > 1 the decay oscillatory. 

In flows subject to CBC, with given vertical and az- 
imuthal wavenumbers, if S > 1 the system supports a 
countably infinite number of inertial normal modes de- 
scribed by Bessel Functions oscillating without growth or 
decay. Otherwise if S < 1 the system allows for only one 
eigenmode - meaning that in general the dynamical re- 
sponse cannot be characterized by discrete normal modes - 
in which case, such linearized disturbances must be treated 
as an initial value problem, a situation well-known in stud- 
ies of plane-Couette flow (Schmid & Henningson, 2000, 
pCf hereafter) . 

We would like to note that the two main results we 
report upon here may be clues in explaining the results of 
Lesur & Longaretti (2005) who show that there is a sub- 
critical transition beyond the Rayleigh line in rpCf when 
SBC conditions are applied whereas no such transition was 
observed for such flows under CBC boundary conditions. 
On a related note, as well, Lithwick (2007) demonstrated 
the existence of steady vortex structures in 3D rpCf flow 
for Keplerian shear profiles when SBC are applied. We 
return to this in the discussion. 



2. Equations, Boundary Conditions and Two 
Integral Statements 

The equations we consider represent the dynamics taking 
place in a "small-shearing box" section located at a cylin- 
drical radius R = i?o centered about the midplane of a 
disk rotating about a central object with the local rota- 
tion vector, fioz. With the external gravity set to zero and 
the density constant, these non-dimensionalized equations 
(Goldreich & Lynden-Bell, 1965,Umurhan & Regev, 2004) 
are 

V • u = 0, (2) 
(dt — qiloxd y )u + u ■ Vu — 2Q,qV — —d x P (3) 
(d t - qQ xd y )v + u • Vm + (2 - q)n u = -d y P, (4) 
(dt — qCloxd y )w + u • Vw = —d z P. (5) 

The flow on these scales is equivalently known as incom- 
pressible rpCf (e.g., Nagata, 1986, Yecko, 2004). In this 
form these equations are the inviscid limit of those consid- 
ered by Longaretti (2002), Mukhopadyay et al. (2005) and 
Yecko (2004). In the language of this paper, x corresponds 
to the radial (shearwise) coordinate, y corresponds to the 
azimuthal (streamwise) coordinate and z corresponds to 
the vertical (spanwise) coordinate - with the correspond- 
ing velocity disturbances u = {u, v,w}. These velocities 
represent deviations over the steady Keplerian flow (as 
manifesting itself in this rotating frame) given by — qfloxy. 
Oo is sometimes also referred to as the Coriolis parameter, 
and in these nondimensionalized units f2o is 1. The local 
shear gradient is defined as the exponent of the general 
rotation law fl(R) = ( J R/i?o) _<z - 

Boundary Conditions. We first present and discuss the 
means by which one administers shearing boundary condi- 
tions (SBC) in the SSB. By this we mean to say that the 
flow variables u, and p are (a) periodic in the azimuthal 
and vertical directions on scales L y = 1 and L z — 1 respec- 
tively and (b) are simply periodic (on scale L x = 1) in the 
shearwise (radial) direction with respect to a coordinate 
frame that moves with the local shear (e.g. Lynden-Bell 
& Goldreich, 1965, Rogallo, 1981). In this work we refer 
to this coordinate frame as the shearing coordinate frame 
(SC for short). The coordinate frame of an observer in 
the rotating frame, will be referred to as the non-shearing 
coordinate frame (NSC for short). In this sense the fun- 
damental equations ([2][5]) have been expressed in the NSC 
frame. 

If (X, Y, Z, T) represent the independent variables in 
the SC frame then we say they are related to the NSC 
variables by X = x, Y = y + qQoxt, Z — z, T = t. If / 
is any dynamical quantity then in the SC frame the peri- 
odic boundary conditions in the shearwise, azimuthal and 
vertical directions respectively are f(X,Y,Z) — f(X + 
L X ,Y,Z), f(X,Y,Z) = f(X,Y + L y ,Z), f(X,Y,Z) = 
f(X,Y,Z + L z ). Of course this means in general that 
f(X, Y, Z) = f(X + L X ,Y + L y ,Z + L z ). As expressed 
in the NSC frame, the periodicity in the azimuthal and 
vertical directions appears like f(x,y,z) = f(x,y + 
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L v ,z), f(x,y,z) = f(x,y,z + L z ), while the periodicity 
in the shearwise direction is f(x,y,z) = f(x + L x ,y + 
qQ tL x , z) . Because of the periodicity in the azimuthal 
direction this statement takes the form f(x, y, z) = f(x + 
L x ,y + y, z) where y(i) — qfl tL x (mod L y ). Note that this 
means that only after exactly integer values of qQotL x / L y , 
(i.e., when y = 0) are the SC and NSC frames coincident, 
meaning every T rop = L y /qQoL x , time units. Evidentally 
the SBC are time dependent when viewed by an observer 
in the NSC frame. 

When enforcing channel boundary conditions (CBC) 
on (PIS]) we require that all quantities are periodic in the 
vertical and azimuthal directions (in the NSC frame) while 
requiring no normal-flow at the channel boundaries. This 
latter statement amounts to u = at x = 0, L x . 

A Pair of Integral Statements. The dynamical equa- 
tions fl2HH]) describe the evolution of velocities which are 
disturbances about a steady rotationally supported flow 
which, on the SSB scales, is the linear shear. Thus, the 
steady solution, u = 0, is interpreted as the undisturbed 
state. We now consider the total velocity v defined as 
v = — qVl^xy + u. As such the governing equations of 
motion ([2][5]) are more concisely written in vector form, 



9 t v + v • Vv = — Vp — 2f2 z x (v + qQ xy), 
V • v = 0. 



(6) 
(7) 



Note that these equations written in this way have ap- 
peared in other works (e.g., Lesur & Longaretti, 2005). 
Taking the inner product of ([6]) with v and following some 
manipulation yields the following 



d t £ + v- V(e+p) = 0, 



— -qttlx 2 . 



(8) 



With use of the incompressibility condition |J7J we may 
integrate ((8]) over the full spatial domain to find, 

^- = - f (e+p)v-MS, £ = [ edV, (9) 
dt Js Jv 

in which V and S are the volume and surface-boundaries 
and n is the unit normal of the bounding surface. £ is 
comprised of (a) the term v 2 /2 which is the kinetic en- 
ergy and (b) the potential-like term — qft^x 2 . The global 
integral £ can change due to the influx of e across the 
boundaries , i.e., f s £~v ■ hdS, and through the external 
work done upon the system denoted by the boundary in- 
tegral f s p~v ■ hdS. The appropriate Bernoulli function for 
this system is e + p. The functional £ is analogous to the 
total energy contained in the flow and its dynamical re- 
sponse (see below) is characterized by the total kinetic 
energy variations in the box. 

By contrast we follow the same procedure as above 
but consider an energy functional comprised only of the 
disturbance velocities. Beginning with (J3]l5j) and taking its 
inner product with the disturbance velocity u, followed 
by the use of ((2|) and either the SBC or CBC boundary 
conditions, one can integrate over the domain to find 



The disturbance energy per unit volume is £. The total 
domain integrated disturbance energy is E = J v £dV. The 
above is known in more general terms as the Reynolds-Orr 
equation (e.g. Schmid & Henningson, 2000) and relates 
the change of the total disturbance energy to the domain 
integrated correlation between u and v. In the absence of 
shear (i.e., q — 0), the disturbance energy is conserved. 

Inspection of the Reynolds-Orr relationship says that 
so long as there is some correlation between the radial and 
azimuthal velocities the disturbance energy must fluctuate 
in time. However from the definitions of £ and E we have 

£ = E+E BheaT ; E shc& , = [ [-qQ xv + \q 2 ^lx 2 } dV, (11) 

Jv 

we see that disturbance energy E is related to £ (the global 
energy) in an intimate way via -Eshear, which we consider 
to be the disturbance energy affected by the shear in a 
dynamical sort of way. 

To appreciate the significance of these relationships we 
first apply CBC to © and find that 



f =0. 

dt 



(12) 



This means that the fluctuations in E come at the expense 
of £? s hear in an equal but opposite way throughout the dy- 
namical response of the flow; i.e., the stresses comprising 
the RHS of (fl7J|) is a measure of the time rate of change 
between these two forms of measured energy. 

We note that if instead of the CBC we require that 
the flow be periodic in y and z (again) but that all flow 
quantities decay sufficiently fast as x —> ±oo, then the 
boundary term in © also vanishes and (fT2")) again follows. 
The same is true if the vertical domain is unbounded and 
disturbances are vertically localized as well, i.e. (fl2|) also 
follows if disturbances similarly decay as z — * ±oo. 

Matters are less clear when SBC are employed because 
£ is no longer conserved, evolving according to ( with 
boundaries of the periodic box set at x = ±|), 



— = qfl [ uvdydz, 
dt Jsi 



(13) 



in which Si/ 2 is the boundary surface at x = 1/2. In gen- 
eral the right-hand side of (fT3]) is not expected to be zero 
and it means that £ can vary with respect to time when 
SBC are employed. The physical reason for this is clear: as 
a fluid parcels exits from one radial boundary with a given 
azimuthal velocity it reenters through the other boundary 
with an azimuthal velocity boosted by the difference in the 
background velocities between the two radial boundaries. 
Another observation about this is that disturbances are 
forced a-priori to have a correlation scale set by the length 
of the box L (Regev & Umurhan, 2008). We note that £ 
will be conserved for SBC if the disturbances are suffi- 
ciently localized in the interior of the computational do- 
main so that there is no power (i.e. fluid motion) near the 
bounding surface Si/ 2 (Nusser, 2007). 
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3. Linear Analysis with CBC 

In considering the response of the flow in the chan- 
nel model we shall assume initially the general form 
f(x, y, z, t) — f af} (x, t) e~sq>(iay + if3z) + c.c. Linearization 
of (J2]{5]) followed by introduction of this spatial Fourier 
form and rearranging results in the single equation for the 
radial velocity 



(d t ~iqn a xa) 2 {d 2 x -f3 2 



a )u a0 — LJ e p u o 



(14) 



in which the epicyclic frequency is defined by lu 2 = 
2(2 — q)flQ. A full breadth of normal-mode solutions of 
the form exp(iLut) exist only for certain combinations of (5 
and a. To see what the restriction is in the following we 
perform a normal analysis and determine precisely under 
what conditions this circumstance occurs. Thus, assuming 
the normal mode form u a/3 — u a/3 exp iut, into (114j) results 



(w - qQ xa) 2 (d 2 - P 2 - a 2 )u a(i + Lu 2 /3 2 u a0 =0. (15) 
The solution to (fT5"|) subject to CBC is given by 

u qj3 = A(lu — qfloxa) 1 / 2 ^ J„ \i\(u) — qCloxa)] Y v [iAaj] — 

Y v [iX(lu — qflQxa)] J„ [i\u>~\ j, (16) 



where 



A = 



ft 2 



2o2 



q 2 n 2 a 2 



v 2 = 



E = 



q 2 Q 2 a 2 ' 



(17) 



where J„ and Y v are Bessel Functions of order v. The 
quantization condition on lu is 

(lu - Woo) 1/2 | J v [i\(tu - Woo)] Y v [iXu] 

-r„[»A(w-Woo)] Ju[iXu]\ = 0. (18) 

The governing equation and its solution has a regular sin- 
gular point for lu = lUqo = qfloa at the point x = 1. 
As it turns out, the solution automatically satisfies its 
requisite boundary condition at x = 1 when lu = Wqo. 
These modes are the fundamental structures that consti- 
tute the so-called wallmodes that give rise to the stratoro- 
tational instability (Yavneh et al,, 2001, Dubrulle et al., 
2005 Umurhan, 2006). 

We observe that if v, as appearing , is a real quantity 
then there is no way for the quantization condition to be 
met by real values of ui (other than the single value Woo) 
according to the known behavior of the Bessel functions 
(Abramowitz & Stegun, 1972). This means that if S < 
1, then there is only one normal mode solution for this 
problem. For S > 1 there exists a countably infinite set 
of normal-mode solutions (overtones) as per the theory of 
Bessel functions (see again Abramowitz & Stegun, 1972). 

Nevertheless, inferring the dynamical content of the 
disturbances, including the shape of the eigcnfunctions 
and the general formulae for the dispersion relationship, is 
not easy. One can more transparently study the behavior 
of these features, including their character change, if an 



asymptotic solution for u afj is developed instead for the 
small limiting values of (3 and a. We proceed by assuming 
£< 1 and, j3 = e(3\, a = eai, followed by the expansion 
procedure lu = euj 1 + e u>s + ■ ■ ■ and u a = u + e 2 u 2 + ■ ■ ■ , 



to the solution of (|T5|) . We find that the lowest order so- 
lutions depend critically on the parameter Si defined by 
Ei = 4:Lu 2 /3 2 /q 2 flQa 2 , which is really a reexpression of E. 
The leading order allowed solutions is uq 



u = A 1 



qQoxai 

LUl 



1/2 



8 In 1 



qfloxai 

LUl 



,(19) 



with 25 = 1/I — Ei and where Aq is an arbitrary ampli- 
tude. The asymptotic form predicts the absence of normal- 
mode solutions if Ei < 1 except for one with lu\ — qftoai 
(the wallmode). The quantization condition on lu\, result- 
ing from imposing the boundary condition at x = 1 results 
in lui = qflo a i/[L ~ ex P( n7r /l^l)], n — ilj ±2, • • •, where 
n labels the particular overtone in question. There are an 
infinite number of radial overtones for any allowable pair 
of /?i's an d ai's. 

In a more general sense: when E < 1, aside from the 
mode lu — lu^, the response of the flow must be consid- 
ered as an initial value problem since the system supports 
no other discrete normal modes. When E > 1 then there 
are two branches of normal modes (i.e., overtones) clus- 
tering about lu = and lu = luoo - in general the fre- 
quencies around which the overtones cluster correspond 
to wavespeeds (i.e., c = Lu/a) that equal the flow speeds 
at the two channel boundaries. 

Modes with n > are the symmetric brethren of those 
with n < 0. The reason is that (fT5|) remains unchanged 
after applying the spatial reflection operation about the 
point x = 1/2 together with the frequency reflection/shift 
lu — > — Lu + qilga. A given eigensolution u^t (x) with eigen- 
frequency lu + generates another eigensolution u^ a J with 



qVL^a such that u a J is the 



mirror reflection of u K ^J about x = 1/2. In general we 



eigenfrequency w_ = — lu^ 

V/3 

see that the approximate theory - including the predicted 
eigenfrequencies - are good even for order 1 values of a and 
(3, but the theory begins to breakdown when they exceed 
order 1 values (Figure 1). 



4. Reflections and Speculation 

1. A Speculation. With respect to the question of dynam- 
ics of wall-bounded rpCf, these results are quite sugges- 
tive. Rincon et al. (2007) showed that Rayleigh-stable an- 
ticyclonic rpCf appears not to show a transition into a 
turbulent state in the manner in which transition is ob- 
served in pCf flow (e.g., Waleffe, 2003). The strategy in 
those problems is to identify steady nonlinear solutions 
of the flow (usually by some continuation method from 
known solutions in otherwise unstable parameter regimes, 
e.g., Nagata, 1990, or forcing methods, Waleffe, 2003) and 
study their stability properties. If there are a collection 
of steady nonlinear structures, some of which unstable, 
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(a) u a p :£ = 19.75, q=3/2 (b) Frequency vs. Horizontal Wavenumber 




Fig. 1. Results from normal mode theory in which E = 19.75 
and q — 3/2 and flo = 1. (a) The exact eigenfunctions for u a p 
depicted for the first three overtones each at two values of a. 
The eigensolutions developed via the asymptotic theory lie ex- 
actly on the eigenfunctions for a = 0.02 where Ao = 1/25. (b) 
The frequencies u> as a function of a for the first three overtones 
are shown (open circles). The predicted frequencies using the 
asymptotic theory are shown with solid lines. The values cor- 
responding to the eigenfunctions depicted in the left panel are 
designated with diamonds. Note that only the positive over- 
tones (i.e., n > 0) are shown here - the negative overtones are 
reflection symmetric about x — 1/2 as discussed in the text. 

then the transition into a turbulent state is thought to 
be approached. These states get easily triggered (bypass 
mechanism) because small scale disturbances in pCf get 
amplified due to the lift-up effect - a mechanism by which 
axisymmetric disturbances result in algebraic instability 
(Ellingsen & Palm, 1975). Rincon et al. (2007) however 
show that there are two problems with this approach when 
employed to identify structures of rpCf near the Rayleigh- 
stable line (e.g., near q — 2 for Re — » oo): that steady 
structures could not be identified and that the anti lift-up 
effect, which is the lift-up effect analog on the Raylcigh 
line, induces radial/ vertical velocity fields that do not un- 
dergo the sorts of secondary instabilities that are charac- 
teristic of the induced azimuthal velocity fields generated 
by the lift-up effect in pCf flow. 

Since the linear analysis of the inviscid Rayleigh-stable 
rpCf shows that for a given a and (3 there are condi- 
tions where there exists a countable number of oscillat- 
ing modes, we openly speculate here that one might con- 
sider identifying oscillating nonlinear structures instead 
of steady ones. Perhaps one strategy would be to (i) iden- 
tify oscillating nonlinear states in the infinite Reynolds 
Number case (Re = oo) and then (ii) ratchet down Re 
and follow/study the structures and their activity until 
either (iii) one identifies a Re of transition, if a unique 
one exists free of hysteretic effects or determine if, as 
Rincon et al. (2007) suggest, rpCf exhibits the charac- 
ter of a chaotic saddle. However we note that identifying 
oscillating (travelling wave) nonlinear solutions is difficult 
in plane Poiseuille flow (pPf) as Toh & Itano (2003) have 



recently shown. However, it might be easier to identify 
such structures in rpCf because of the presence there of 
linear oscillations whereas such oscillations are absent in 
linearized pPf. 

2. Some Reflections. An individual wavemode subject 
to SBC, i.e. ~ exp (ikX + iaY + i(3Z), where k, a and 
(3 are wavenumbers as before, was shown by Balbus & 
Hawley (2007) to be an exact solution of the nonlinear 
equations of motion (i.e. [2][5]). This is because this SBC 
waveform ansatz, together with the incompressibility con- 
dition, results in the nonlinear terms exactly cancelling to 
zero for each velocity component. Given () 1 Ofl 1 1|) and (|13jl 
it follows that for these solutions the following holds, 

d£ _ dE 
~dt ~ ~dt' 

where E is the disturbance kinetic energy of the wavemode 
(see ITOf . From (|11[) it follows that these disturbances con- 
serve the shear energy ©shear- A study of the form of .©shear 
shows that this is a natural outcome due to the symme- 
try of the modes with a,/3 0. However, for modes in 
which a — /3 — 0, an inspection of the governing equa- 
tions of motion in the SC shows that it necessarily fol- 
lows that the corresponding azimuthal velocity mode, i.e. 
which is uniform in the vertical and azimuthal directions, 
is zero. Thus, by its very construction, these exact non- 
linear modes as a result of SBC cannot support nonlin- 
ear modifications to the shear. Consequently, such alter- 
ations can occur only via three-wave interactions (Balbus 
& Hawley, 2007). We note that no such disturbances are 
supported in CBC flows - there always exists some amount 
of self-interaction for most modes of the system. 

We note that for a general collection of perturbations 
utilizing SBC, (fl3|) says that £ ^ 0. The fluctuations in the 
total energy relate directly to the lateral exchange of ki- 
netic energy from radially neighboring copies of the shear- 
ing box, where azimuthal velocities enter or exit boosted 
by the difference in the background shear velocities be- 
tween the two boundaries - the quality of this exchange is 
a result of the type of correlations imposed upon the dis- 
turbances. For SBC it means that these correlations are 
set a priori by the periodic length scale of disturbances 
which is set to be on the scale of the box L (e.g. see Regev 
& Umurhan, 2008). By contrast, the imposition of no- 
normal flow boundary conditions removes this source of 
energy into and out of the domain for problems invoking 
CBC, thus £ = in that case. 

We feel that these observations about the differences 
between CBC and SBC, together with recent results, 
raises more questions that need to be resolved. For in- 
stance, linear disturbances of the CBC type are under- 
stood not to result in decay because of reflections off of 
the channel walls. SBC conditions, on the other hand, sup- 
port modes which can result in fluctuations of the total 
energy depending on the correlations that exist across the 
radial boundaries. It is interesting, therefore, that Lesur 
& Longaretti (2005) demonstrate the persistence of activ- 
ity for simulations of the SSB (slightly below the Rayleigh 
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line) using SBC while they report the presence of no such 
activity for CBC. On the other hand, there is in Shen et 
al. (2006) a recent demonstration of the decaying quality 
of disturbances in the SBC for flows run with q = 3/2. 
Equally intriguing is that Lithwick (2007) shows that 
a single vortex may persist for a long time in a three- 
dimensional simulation of the SSB for q = 3/2 using SBC 
which involves the complex interplay of three-wave mode 
interactions. It is interesting to us that while CBC con- 
serves its total energy £ no activity has yet been observed 
under Keplerian-like conditions while for SBC, which does 
not conserve its £ , sustained activity occurs in some cases 
while decay characterizes others. 

Given the preceeding discussion it seems there are im- 
portant questions to be asked: 

1. What is the role of the non-self interacting wave modes 
resulting from the use of the SBC in the maintenance 
and/or destruction of activity in SSB simulations? 

2. What influence does the prescribed velocity correla- 
tions, which result in the total energy not being con- 
served (cf. (fl3|) ) have on the overall evolution of general 
SSB simulations using SBC? 

3. What influence does the non-self-interacting nature of 
these waves in the SBC have on the velocities on the 
radial boundaries of the system which, in turn, affects 
the amount of energy entering/leaving the system.? 

4. How are we to understand these effects in relation to 
dynamics supported in the nonlinear response of flows 
utilizing CBC (i.e. wall-bounded rpCf)? 

These observations and questions underscore the need for 
further analysis (also see Dubrulle, et at, 2005). 
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